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Abstract 

We consider the one dimensional focusing (cubic) Nonlinear Schrodinger equation 
(NLS) in the semiclassical limit with a one parameter family of exponentially decaying 
in absolute value initial conditions. We prove smooth parametric dependence of the 
asymptotic solution on the parameter utilizing the Riemann-Hilbert approach. Nu- 
merical computations supporting our estimates of important quantities are presented. 

1 Introduction 

We study a type of scalar Riemann-Hilbert problem that appears in semiclassical or small 
dispersion limit of integrable systems. We prove the smooth dependence of the solution on 
a crucial parameter. Our motivation comes from the semiclassical focusing NLS equation 

la ma mi us] 

ied t q + e 2 d 2 x q + 2\q\ 2 q = (1) 

with the initial condition 

q(x,0) = A(x)e^ s{x \ A(x) > 0, /i > (2) 

in the limit as e — > 0, where the parameter /i is the focus of this paper. For any value of 
e > 0, the solution process requires solving a 2 x 2 matrix Riemann-Hilbert problem (RHP) 
in the complex plane of the spectral parameter z of an underlying Lax pair operator [231 EI] • 
The quantities x, t, and /i enter as parameters. The Riemann-Hilbert approach is used 
to other integrable systems in general as established by [19] and it is a major tool in the 
asymptotic analysis of integrable systems as established with the discovery of the steepest 
descent method [TQ1 E]. The asymptotic methods via the RHP approach also apply to 
orthogonal polynomial asymptotics [HI [7], and to random matrices [U EJ [121 113] . 

Riemann-Hilbert problems formalism is another form of the (2D) potential theory in the 
complex plane of the spectral parameter z. The main question is to find the equilibrium 
measure of a set - a minimizer of an energy functional [61 [20] . This equilibrium measure is the 
spectral measure of the underlying linear operator. In the case of NLS, it is (non-selfadjoint) 
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Zakharov-Shabat system (a Dirac-type operator). For self-adjoint operators the equilibrium 
measure is supported on a subset of the real axis. In a general non- self-adjoint case the 
support is a subset of the complex plane. A typical support of the equilibrium measure in 
the case of NLS in the semiclassical limit seems to be a finite union of arcs in the complex 
plane with complex conjugate symmetry jH [Sj, [16], |23j [25] . Denote the end points of the 
arcs as {otj}j =0 with some finite N' G N. This leads to finite genus Riemann surfaces in the 
asymptotic analysis. 

There are several approaches to obtaining asymptotic solution from the initial data. 
One approach is to analyze the RHP for a fixed x and t. Another approach is to look at 
the dependence of the RHP on x and t as parameters. This includes the analysis of the 
contributing limiting arcs through analysis of the end points ocj. These are often referred as 
branchpoints since they arise as the branchpoints of a square root. 

The smooth dependence of the branchpoints aj with respect to parameters x and t was 
obtained in [21] through formulae of the form 

2m^(a j ,a,x,t) 
dx [X ' t] ~ D (a,x, t) f rt ^ff* <*C ' [ ) 

\ ' ' / J7 (C-aj(x,t))R(C,a) >> 

do^ = 2m^(a j ,a,x,t) 
dt {,! Did x t)S /,(c,a! ' f) dC 

We consider initial data being a semiclassical approximation [23] of 

g(x,0) = -sech (x)e-^^ tanh(s)ds , (jt > 0. (5) 

This family of initial conditions is interesting because of a transition at fi = 2 with solitonless 
interval (fi > 2) and radiation with solitons interval < fi < 2. It has been studied in a 
number of papers [2"3"| |2"H |2"B] . The semiclassical limit has been completely analyzed for \x > 2 
while for the soliton case < fi < 2 the answer is known only for some finite times. The RH 
approach runs into difficulties with the error estimates and was not able to continue past 
a certain curve. Numerical experiments have shown absence of any noticeable transition in 
the behavior of the arcs end points at the critical value // = 2 [3]. The perturbation 



theorem for NLS 3.6 establishes this fact rigorously. 

We extend formulae (|3]|4]) to the case of the contour explicitly depending on an external 
parameter /1, 7 = ^(fi) 

[x, t, n) - /'(C,*,*,m) dC {b) 

and the dependence is smooth, meaning that the contour, the jump matrix, and the solution 
of the scalar RHP evolve smoothly in \x. Moreover we simplify the expression for |f as J37b. 



We also compare (see Fi g. [3] ) formula (|6| with the direct computations of the branchpoints 



from solving the system (13). 



Moreover, we prove the preservation of genus of the solution for all \i in an open interval 
around any fj, Q for which the genus is known. In particular, the genus is preserved (0 or 2) 
for all x and t > for some open interval (which depends on x and t) for \i < 2. 
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Figure 1: The RHP jump contour in the case of genus 4 with complex conjugate symmetry 
in the notation of |23|. 



The paper is organized as the following: Section 2 - main definitions and prior results 
are stated, Section 3.1 - we prove analyticity of / in /i, which leads to differentiability of the 
branchpoints <x,- = aij((j), Section 3.2 - we consider \i dependence of all main quantities stated 
in theorem 3.6, Section 3.3 - we additionally consider sign conditions and prove preservation 



of genus in theorem 3.15 



2 Preliminaries 

We consider the scalar Riemann-Hilbert problem (RHP) which arises as a model (simplified) 
RHP in the process of asymptotic solution of the semiclassical focusing NLS ([TJ) with the 
initial condition (|5| 

h+(z) + h_(z) = 2Wj, on 7 m>i , j = 0, 1, N, 
h+(z) - h-(z) = 2%, on j_ cd , j = 1, ...,N, (7) 
h(z) + f(z) is analytic in C\7, 

where the jump contour 7 belongs to the class of admissible contours with the following 
properties. 

Definition 2.1. 

We say 

7 G r(a,/z) 

if a simple contour 7 consists of a finite union of oriented finite length piecewise continuous 
simple arcs in the complex plane 7 = (U7 mj ) U (U7 CJ ) with the distinct arcs end points 
a = {ctj} 4 ^ 1 , as shown in Fig. 1. We denote the main 7 mj - and the complementary arcs 
7 CJ - as the following 

7m,o = [ai, Oia] , J m ,j = [«4j-2, a 4 j] U [ocaj+i, «4j-i] , j — 1> •••) -^V> 
7c,j = [«4i-4, a 4 j_ 2 ] U [a 4 j-i, aty-3] , j — 1) •••) N. 
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(a) Contour 7 (b) Contours 7 



Figure 2: Contours of integration for function h(z) (J8J). Zo is a point of non-analyticity of 
f(z) on 7, and = | depends of an external parameter /x. 



Let zo = 7} G 7, 6e a point 0/ logarithmic singularity of f. 

For simplicity we assume that 7 is not tangent to the real axis at z = | > 0. We 
additionally assume that the contour 7 has complex conjugate symmetry 7 = 7. 

Note that for fixed a and /x the contour 7 aside from passing through z = ctj and z = | 
is free to deform continuously within of domain of analyticity of /. Thus for a fixed /, a 
contour 7 G r(a, /x) depends on the arcs end points a and on /x through z , 7 = 7(0, xx). 

An immediate fact about these contours is the following. 

Lemma 2.2. Let 70 = 7o(«o,A i o) G r(a ,/Xo) wii/i the distinct arcs end points a . Then 
there exist open neighborhoods of c? and /io such that for all a in the neighborhood of Oq, 
for all fi in the neighborhood of fio there is a contour 7 (a, /i) G r(a, /i). 

Definition 2.3. 

VKe say 

7 G f (7,5,/i) 

if 'y is a simple closed contour around 7 G r(a, /i) within the domain of analyticity of f . 
with complex conjugate symmetry 7 = 7. 

Remark 2.4. Similarly we define j m j and 7^-. 

Remark 2.5. By considering the loop contours 7, 7 m ,j, 7cj, the explicit dependence of the 



contours on the end points a is removed (for example in (29 32 )). So even though 7 = 7 (a, //), 
in all our evaluations below 7 = 

Remark 2.6. Lemma 



2.2 



implies that if 70 G f (70, a , /xq) then there is a contour 7 G r(a, ji) 



such that 7 G T(7,a, /i) for all a and /x in some open neighborhoods of a and /x - 
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The solution of the RHP (J7|, h(z) can be found explicitly 



h(z) 



2m 



j=0 J ~lm,j 



(C - z)R(C) 



N 



1 •'Tcj 



(C - *)*(0 



or in the determinant form 



h(z) 



D 



K(z), 



where z lies inside of 7 and outside all j Cj j and jmj, and where 

J7m,l #(C) ' ' ' •'7m,l -R(C) J7m,l 



and 
with 



1 

2m 



J*/m,n R(Q 

£ 



r _d(_ 

J%,n R(0 

/(CK 



R(0 
£ C^C 

J 7c,i 



7m, n R(0 

N ~ 1 dC 
R(0 

7c,n i?(C) 



.4 



F 7 R(0 




h R(0 




D = 


det(A), 


( f 

J 7m,l 


d( 
R(0 


£ i N ~ ld < 
' ' ' Jj m ,i R(C) 


J 7m.n 

J 7c,i 


d( 

R(0 
dC 

R(0 


JX lm,n R(Q 

' ' ' hc,i R(0 




dC 
R(0 


£ C N ^dC 

' ' ' J%,n R(C) 



(C-*)R(0 
£ rfC 

Jjm.n (C-Z)R(C) 

£ d $ 

r%,i (c-z)r(o 

£ rfC 
hcN (C-z)R(O 

h (c-*)fl(C) 



The arcs end points {otj} satisfy the system 



(9) 



(10) 



111 



(12) 



K(a n 



0, J = 0,1, 



AN + 1. 



(13) 



The dependence on x and t was considered in [21]. This is a simpler situation when the 
jump contour 7 in the RHP §7§ is independent of the parameters. 
The main related results in [21] are the determinant form ^ and 

Theorem 2.7. Let f(z,(3), where (3 E B C IR m . For all (3 E B assume f(z,(3) be analytic 
in on S E C. Moreover, j\S consists of no more than finitely points and f continuous on 
7. The modulation equations JTgp imply the system of AN + 2 differential equations 



(a 



no 



7 (C-aj)fl(C) 



(14) 
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In particular, one gets ^ and Q for parameters x and t. Note, that the contour 7 is 
assumed independent of parameters x and t explicitly. The dependence on these parameters 
comes in through the branchpoints a = a(x,t). 

The main related result in [22] is 

Theorem 2.8. Let the nonlinear steepest descent asymptotics for solution q(x,t,e) of the 
NLS |7]j be valid at some point (x ? to) ■ If (#*,£*) is an arbitrary point, connected with (x , to) 
by a piecewise- smooth path S, if the contour 7(2, t) of the RHP ([?]) does not interact with 
singularities off(z) as (x,t) varies from (x ,t ) to (x»,t*) along S, and if all the branchpoints 
are bounded and stay away from the real axis, then the nonlinear steepest descent asymptotics 
(with the proper choice of the genus) is also valid at (x*,t*). 
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We extend Theorem 2.7 and make partial progress in the direction of Theorem 2.8 
the case when the jump contour explicitly depends on the parameter \x. We require that 
the point of logarithmic singularity of / = f(z,fi), Zo = | is always on 7. Additionally 
we prove preservation of genus for all x > 0, t > 0, /i > 0, under certain conditions which 



guarantee that the parameters are away from the asymptotic solution break Theorem 3.15 
In particular, the genus is preserved in the neighborhood of the special value of the parameter 
fi = 2. Thus we obtain that for all x > 0, t > (except on the first breaking curve) there is 
a small neighborhood such that for all fi < 2 in the neighborhood, the genus is the same as 
for jj, — 2, where it is known to be or 2. 



3 /i-dependence in the semiclassical focusing NLS 



3.1 Setup 

To apply the methods from |21j we need analyticity of f(z,n) in the parameter /1. 

The function f(z) comes from a semiclassical approximation of the family of initial con- 
ditions for NLS Q as [23]: 



f(z,fi,x,t) 



— z 



™+ln(^ 
2 V2 



+ In (z + T) + —— In {z - T) 



Ttanh 



-1 



T 
/V2 



xz — 2tz H — In 2, when Qz > 
2 



and 



f{z) = f(z), when Qz < 0, 



(15) 



(16) 



where the branchcuts are chosen as the following: for < fi < 2 the logarithmic branchcut 
is from ^ along the real axis to +00, from T to and along the real axis to +00, from —T 
to and along the real axis to —00; for \i > 2 - from T to +00 and from —T to —00 along 
the real axis, where 



T = T(/i) 



El. 
4 



ST > 0. 



(17) 



For > 2, T > is real and for < fi < 2, T is purely imaginary with ST > 0. 



6 



Then for ^sz > 



Z,fj) 



iri 1 rii 
— + - In ' 



z) +ln2 + 



8T 



hi(z + T) — ln(z — T) — 2 tanh' 



.1 2T 



where tanh 1 x = x + 0(x 3 ), x — > then for 7^ 



<9/i 



^ = J + \ ln ( f " z ) + ln 2 + Tz ~ 5 + 0(T) ' T 



1 



So /i = 2 is a removable singularity for fn(z,fi) and 

C? f . , 7T2 1 , . , 1 1 

lim ■ 7 ^-(z,u) = h - n(l - z) +ln2H . 

M W^ 1 ,W 4 2 v ; 2z 2 



(19) 



(20) 



which is analytic in z for Qz 7^ 0. 

Lemma 3.1. f(z, ji) and f'(z, /i) are analytic in fi for \i > 0, x > 0, t > 0, for all z, Qz 7^ 0, 
z£[-T,T\. 



Proof. Consider 



Til 



hi 



z ) + - In [ z 



2 V 



+ 1 - x - Atz, 



(21) 



which analytic in \i > 0, for ^ 0, z ^ [—T,T]. 

For > 0, n 7^ 2, f(z,/i) is clearly analytic in /1 for 3z 7^ 0. At /i = 2 (T = 0) we find 
the power series of f(z,fj,) in T and show that it contains only even powers. Since 



-,2k ( A 4 



- 1 



it will show analyticity of f(z,fi) in \i. 
Start with expanding in series at T = 



0* + 2)' 

4fc 



(A* " 2) i 



1 00 00 

^vW^c^, ln(,±T) = ^ 



k=0 



n=l 



n 



Then the terms in (15) are analytic in // 



2 — T 

ln (z + T) + ln (2 - T) 



--e K§)"-£J(?)" 



n «s etien 
00 

'In 2 + 



' 2Ac(2fc - l)z 



2&-1 



T 



2k 



(22) 



(23) 



(24) 
(25) 

(26) 
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which has only even powers of T and is analytic in /x for ^sz ^ 0. The inverse hyperbolic 
tangent term in (15) and taking into account that tanh - z is an odd function 



rri rri C- 1 

VManlr 1 — = Ttanlr 1 , = Ttanh" 1 ^ c fc T 2fe+1 



M/2 



vTTt2 



k=0 



n2k+2 



k=0 k=0 

which also has only even powers of T. 

So f(z,x,t,/J l ) is analytic in /x for /x > 0, x > 0, t > 0, Qz ^ 0, z [— T(/x), T(/x)]. 



(27) 
(28) 

□ 



Note the jump of f(z) is caused by the Schwarz reflection (16) on the real axis and it is 



linear in z since 3/ is a linear function on the real axis (as a limit) near | with Qf f |J = 

ESI. 



3.2 Parametric dependence of scalar RHP 

The main difficulty is the dependence of f(z) (thus the RHP 0) and the modulation equa- 



tions (13) on parameter /x which also controls the logarithmic branchpoint z = | on the 



contour 7. We show that the dependence on /x is smooth. 

To solve K(a,n) = 0, we need the Implicit function theorem and nondegeneracy of the 



system. The following technical lemma simplifies the expressions (j8j), (10). 
Lemma 3.2. 



Let function f be given by (15) and there is a contour 7 = 7(0, /x ) E T(a, /x ) has fixed 



arcs end points a. Then there is an open neighborhood of /xo such that for all /x in a small 
open neighborhood of Hq there is 7(/x) E T(7,a, /x) 



d_f C7(C^R 



c 



d 



I 



a 



7O) 



t(m) 



i?(C,a) 
(C-aiWC,5) : 



n G N, 



n G N, 



(29) 



(30) 



d_ f C n d( 







a: 



0, j = l,2,...,xV, nEN, (31) 



9// j^ cj (C - &j)R((, 



a) 



0, j = l,2,...,N, nEN. 



(32) 
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Proof. 



The idea of the proof is to consider finite differences and take the limit as A/i — > 0. The 
complication is that both the integrands and the contours of integration depend on /a. 



Denote the integral on the left in (29) as I± 



w-L^*- (33) 

where j(fi) G r(7, a, /a). Consider 

h{fi + A/i) -Ji(/i) 
A/i 

with small real A/i 7^ 0. There are two logarithmic branchcuts near the contours of integra- 
tion: in f(z,fi) and in f(z,fi + A/i) with both branchcuts are chosen from z (fi) := | and 
zq{h + A/x) horizontally to the right along the real axis. Additionally, these functions have 
a jump on the real axis for z < ^ from Schwarz symmetry. 

We choose some fixed points Si and 62 to be real, 5i < f — < f < Both 

contours of integration 7(/i), 7(/i + A/i) are pushed to the real axis near zq and split into 
[Si, S 2 ] '■— [Si + zO, <5 2 + iO] |J [5 2 — zO, <5i — iO] and its complement. On the complement, we 
can also deform both contours to coincide. So 7(/i + A/i) = 7(/i). 

Note that across [^i,^], /(z, /i) has a jump 7u |zo(a0 — z\ and f(z,fM + A/i) has a jump 
7rz |zo(/i + A/i) — z|. So contributions near zq in both cases are small. 

Then 

hiji + Ari-hQi) _ 1 (i Cf((,» + A/i) rfC _ r Cj%A dC \ (34) 



A/i A/i \J^ +A ^ R{(,a) J w R{C,a 

we add and subtract <£, > - ^/^t^- dC 

1 ^ c"/(c,/. + m (JC+ / c"(/(c,M+Art-/(c,A.)) dc| (35) 



A/i V/7(h+Ah)-7(m) R-iCl®) />(«) -R(C)«) 



The first integral is 0, because 7(/i + A/i) = 7(/i). 
Thus 

A// R((,a) 1 ; 

The last step is to take the limit as A/i — > and to interchange it with the integral. The 
contour of integration is split into two: a small neighborhood near zq and its complement. 
For the integral near zq, by a direct computation can be shown that the limit can be passed 
under the integral. The integral over the second part of the contour has the integrand 
uniformly bounded in \x since log — ^) in |^ is uniformly bounded away from ^, so the 
limit and the integral can be interchanged. This completes the proof for the first integral 



(29). 



The second integral (30) is done similarly. The rest of the integrals (31)-(32) are inde- 
pendent of /i since the only dependence on \i sits in z (fi) G 7 m ,o- 

□ 
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Using lemma 3.2 









• f 








"7m, 1 R(0 


J 7ra,l 








£ 






dK _ 

— [a^a.y) = 


1 

2m 


r d( 
J%,i R(0 ■ ■ 

£ 


J 7c,i 


R(0 






J%,n R(0 " " 










h R(0 


■ *<" 


_1 /m(CK 








fl(0 



J7m,l 
J 7m,n 

J 7c,i 



(C-«j)«(C) 

(C- aj )R(0 
(C-«j)fl(C) 

7civ (C-oy)ii(O 
J't (C-aj)Jl(C) 



(37) 



where / M is given by (18). 
Lemma 3.3. 

Let f be given by (15) and the contour 7 = j(do,fj, ) G r(a ,/xo) has the arcs end points a . 
Then there are open neighborhoods of a® and /x such that for all a and /x in the neighborhoods 
of ao and /xq respectively there is a contour^ = j(a, /x) G T(a,fi) and 



j = 0,l,...,4iV + l 



(38) 



Kj(a,n) := K(a,j , a, /x) , 
is continuously differentiable in a and in /x. 
Proof. 

First we notice that since 7 G r(a 0) A t o) by Lemma 2.2, there are neighborhoods of a and 
Ho such that for all a and /x in the neighborhoods of do and /xo respectively there is a contour 
7 = 7(«>^) e r(a,/i). 

Kj(d, /x) is analytic in a by the determinant structure and the integral entries ( 10 ), where 
explicit dependence on a is only in the R(z, a) term which is analytic away from z = a 

j (' i v • 1- 



31 



The integrals in the last row of fr^oy, a, /x) in (37) involve the function f^ given by (18) 
which is integrable near z — | and hence fr^ay, c?, /x) is continuous in /x. Thus Kj(aj, a, /x) 
is continuously differentiable in a and in /x. 

□ 



By Lemma 3.3 the modulation equations ( | 1 3| ) 

Kj(d, /x) = K(aj, a, xx) = 



(39) 



are smooth in a and in the parameter /x. Next we want to solve this system for a = a(/x) 
and conclude smoothness in /x by the Implicit function theorem. 



For the next lemma we need K'(z, a, /x) = ^r{z, a, /x) 



K'(z, a, /x) 



2ttx 



J7m,l «(0 



fid BlU 



7m, iV -R(C) 



7c,;v H(C) 
/(CK 
7 fl(C) 



-„..j «(C) 

'7m,iV R(0 



^7m,l (C-^) 2 «(C) 
■^7m,iV (C-^) 2 ^(C) 

£ rfc 
J 7c ,! (C- Z ) 2 i?(C) 



7c,iv fl(C) 



7c,jv (C-^) 2 «(C) 
x /(CR 
3*7 (C-*) 2 -R(C) 



(40) 
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where z is inside of *f(fj) and inside of j m j and %j or 



Lemma 3.4. 

Let f be given by (15) and a contour 7 G r(c? ,/xo); where a and /x satisfy 

K(a ,no) = 0. 

Assume that for a = {o^\ , lim^^o K'(z, a , /x ) ^ 0, j — 0, 1, . . . , 4iV + 1. 
Then the modulation equations 

K(a,fj.) = 

can be uniquely solved for a = c?(/x) which is continuously differentiable for all /x in some 
open neighborhood of /iq and a(/x ) = So- 



Proof. K is continuously differentiable in a and in /x by Lemma 



As it was shown in [21], the matrix 



dKjoj 
dctj 



dK 



h(z) 



3,1 



dK{aj) 
dan 



3.3 



is diagonal and 



So 



-D lim . 

2 *->«o V.i?(z) 



lim a, /x) ^ 0. 



det 



da 



n 



<9a,- 



^0 



(41) 



(42) 



under the assumptions. By the Implicit function theorem, d?(/x) are uniquely defined in some 
neighborhood of /x and smooth in /x. Note a(/Xo) = «o by assumption. □ 

Remark 3.5. 



The condition lim^^o ^'(z, c?o, /^o) 7^ 0, j = 0, 1, . . . , 4A" + 1 in Lemma 3.4 is equivalent to 
0,^ = 0, l,...,4iV + l. 

All quantities below depend on parameters x and t. We assume that for the rest of the 
paper x and t are fixed. 

Theorem 3.6. (/i-perturbation in genus N) 

Consider a simple contour 7 = 7(/xo) consisting of a finite union of oriented simple arcs 
70 = (LWmj) U (U7cj) with the distinct arcs end points do and depending on parameter /x 
(see Figure [l]). Assume do satisfies a system of equations 

K (a , /x ) = 0, 



and f is given by (15). Let 7 = 7(0, /x) be the contour of a RH problem which seeks a 
function h(z) which satisfies the following conditions 



h+(z) + h-(z) = 2W j , onj m ,j, j =0,1,..., AT, 
h+(z)—h-(z) = 2Sl j , onjcj, j = l,...,N, 
h(z) + f(z) is analytic in C\j, 



(43) 
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where Qj = Qj(a,fi) and Wj = Wj(a,n) are real constants whose numerical values will 
be determined from the RH conditions. Assume that there is a function h(z, c?o, /zq) which 



satisfies (43) and suppose fc j^'"^ - 1 7^ for all z on 7. 
Then the solution a = a(fi) of the system 



K (a, fi) = 



(44) 



and h(z,a(n), n) which solves (43) are uniquely defined and smooth in n in some neighbor- 
hood Of jlQ . 

Moreover, = Qj(a(fi), n) , and Wj(^) = Wj(a([/,), /x) are defined and smooth in jj, 

in some neighborhood of fiQ. 
Furthermore 



dacj 
dfx 

dh 
djj, 



(aO 



z,d,ji) 



2m 



dfj, 



D{a{ii),n)§, 
R(z, a 



no 



where z is inside ofj, 



2m 



7 ( M ) (C-o,fci))Jl(C,30*)) 



%,aO 



7(M) 



(C-z)R(Ca 



(45) 
(46) 



9Wi 
djj, 



D 



D 



£ ^ 

Jjm,l R(0 


Jy m ,i R(0 


£ JZ- 

JX / m ,N R(0 

£ ^ 

J%,i R(0 


£ C N -'dC 

- r 7m,iV R(0 

£ ' C^-^C 
J7cl fl(C) 


£ ^ 

J%j-i R(0 ' 
J-y R[C,a) * 

hcj+i R(0 ' 


x c^-^c 

• J 7cj -i R(0 

■ J7 R{(,d) a >> 
£ C N ~ ld C 
' R(0 


®7c,iv K(C) ' 


£ ( N ~ ld C 
J%,n R{0 


£ d $ 

J7m,l H(C) 


X C*" 1 * 

J*tm,l R(0 


£ 

J7 mj _i R(C) 

J 7mj + 1 -R(C) 


x c"- 1 ^ 

J7mj-1 -R(C) 

£ < W>dC 

■ ■ J 7 R(CA) ^ 

x c"- 1 * 

J 7mj+1 R(0 


£ J*L 
Jt'to.jv -R(C) 

J%l H(C) 


x c^-^c 
x c^-^c 


r dc 


x t N -;d( 

J%,n R(0 



(47) 



(48) 



where R(C) = R((,dfy- 
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Proof. By Lemma 3.4 
similarly as in [21 



are continuously differential) le in /i. Formula for ^ is derived 



which define a = a(/x) with respect to fi 

AN+l 



We differentiate the modulation equations K(aij) = K(aj,a, /x) = 
t to fl 

dKiptj) dai OK 



dai djJL <9/x ^ aj ' 



(49) 



where the matrix j 8J ^° J - 1 j is diagonal [2TJ so 

dK(a.j) dctj dK(aj] 



Since 



dK(aj) D(a,fx) 



we arrive to the evolution equations for af 



2 ™f « 



d» D{a,n)§ f 



(C- aj )R{C,S) 



j = 0,...,AN + l. 



Next we compute p-(z,a,fi) which satisfies the scalar RHP 



Then 



9n,+( z ) ~ 9»,-{ z ) = U( z )i z e 7m j, j = 0, 1, N. 



V*' a,AiJ ■ 2™ L((-z)R((,a) 



(50) 
(51) 

(52) 

(53) 
(54) 



where 2 is outside of 7, §^(2) behaves like log(z — z ) near z = z = |. So a, /x) satisfies 



(46) where 2 is inside of 7. 

Constants Wj and Qj are found from the linear system [21 

" - ~ N - ( n W . 



Differentiating in /x and using Lemma 3^ leads to 



•d£ = 0,n = 0,..JV-l. (55) 



C7,(C,aO 



/ 7 ( M ) a 
or in matrix form 

l "W K(C) 

r _dC_ 
J7cl i?(C) 

r _d(_ 



•X7m,;v -R(C) 
J 7c,i 



/ avy \ 




[ dfi u 




\ an 




\ dp k J 





( i 



jv-i 



(57) 



/ 
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So |^ and ^f- satisfy (47) and (48). Note that D ^ for distinct a/s [2T 



□ 

Remark 3.7. In [21] was considered the case when the contour 7 was independent of external 
parameters x and t. We consider a case of the dependence on /i when the jump contour 
explicitly passes through point of singularity of /. 



Remark 3.8. The perturbation theorem 3.6 guarantees that the solution of the RHP (43) is 
uniquely continued with respect to external parameters. Additional sign conditions on ^sh 
need to be satisfied, for h to correspond to an asymptotic solution of NLS as in [23]. The 
sign conditions have to be satisfied near 7 and additionally on a semiinfinite complementary 
arcs connecting the arcs end points of 7 to 00. 

3.3 Sign conditions and preservation of genus 

The error estimates in the RH approach to semiclassical NLS require certain sign conditions 
to be satisfies. In particular ^sh(z) = on 7 m j and ^sh(z) > on j C) j. 
On the real axis (in the limit from the upper half plane), 

3f( z + iO) = lim f(z + i6), zeR. 
is a piecewise linear function [23] 

9/<«+«)-{jjf:g!- ;>\; 

and since g(z) is real on the real axis, Qh(z + iO) = —Qf(z + iO). It is important for us 
that the growth of ^sh(z) can be estimated and be bounded from below away from zero as 
z — > oo. 
Similarly, 

( f ^<o, 
9f/'(z + i0) = I -f < 2 < f, 

I 2 Z > 2' 

and since g^- 2 ) i s rea l on the real axis, Qh'(z + zO) = —Qf'(z + iO). 

Definition 3.9. Define 7 00 = 7°°(a, //) as an extension of a contour j(a, /i) G /i) 
as 7°°(a, / u) = (00, a^] [J 7(6!, /x) U^AM-i) °°)- additional arcs are considered as 

complementary arcs 7 Cj0 = (00, c^n] C C_, 7 c ,jv+i = (00, c^at] U[ a 4AM-i5 °°) an< ^ assume 

7c,N+l = lc,N+X, SO 7 00 = 7°°. 



Lemma 3.10. If £ne conditions of Theorem 3.6 holds on 7°°(6!o,/io) for K(ao,Ho) = then 



the statement of the theorem holds on j°°(d,fi), where K(a,fj,) = 0. 
Proof. 

The proof is unchanged since / is analytic near the additional semi-infinite arcs 7 Cj o and 
7 Ci at + i and the jump condition for on the additional arcs 7 Cj o and 7 c ,at+i is taken to be zero 
(do = 0, tt N+1 = 0) [23]. □ 
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Note that the conditions in Theorem 3.10 are more restrictive since 7 C 7 



Definition 3.11. A function h satisfies sign conditions on 7 00 if^sh(z) = if z G j m> j and 
S*h(z) > if z £ j Ctj for allj = 0,...,N+1. Denote h G SC(y°°). 

Note that the sign conditions on 7 mj - are satisfied automatically (^sh(z) = 0) through the 



construction of h(z) by (J8J) in the case of h solving a RHP (43). We only need to check the 
nonnegativity of on the complementary arcs, especially on the semi-infinite arcs (00, a 4 jv] 
and [ck4jv+i, 00). 



Theorem 3.12. Let f be defined by (15). Let K(a ,jj, ) = 0, 7^ G T(ct ,jj, ) and assume h 
solves RHP(^, a , // , /) with y^ggy? for all z G 7^°, and h G SC{y™). 

Then there is an open neighborhood of fi^ where for all fi, there is an h which solves 
RHP(j°°,a,/i,f) with 7 00 = 7°°(a,/x), K{a,n) = 0, ^ for all z G 7°°, and h G 



Proof. 



Take any /i in a small eno ugh o pen neighborhood of /zq- There are two things we need 



to prove in addition to lemma 3.10 the sign conditions on 7 00 and that ^ on 7 . 

Assume h G SC (7^°). Then h G SC^di)) by continuity of h in z and /1 and compact- 
ness of 7. The semi-infinite arcs (00, a^} and [a4jv+i, 00) we can pushed to the real axis as 
(-00 + iO, -f + iO) |J [-f + iO,a 4N ] and (-00 - iO, -f - iO) [j [-| - z'0, a 4JV +i] respec- 
tively. On (— 00 — 2O, — ^ — z'0) and (— 00 + iO, — | + z'0) , Qh(z) is positive and [— ^,a 4 Ar] 
and [— #,Q!4jv+i] are compact. So > on 7°°(/i), that is h G S , C(7°°(/i)). 



Let ^ ( - z .'. Mo \. ^ on 7^°. Then there is a constant C > such that 



h'(z,n ) 



R{z,a(ti )) 



> C for all 



2 G 7^°. Consider the solution h(z,/j,) of RHP^ 00 , a, where K(a,n) = 0. By theorem 
3.6| and lemma 3.10 such function exists and continuously differentiable in /1. Moreover, 



h'(z, jj) is continuous in /1. Since 7 is a compact set in C and R r^^vj is continuous in z and 

fi, we have ^g^) ^ for a11 z G 7 ' 

Now as before, the two semi-infinite pieces (—00 + iO^^] and [«4jv+i, — 00 — iO) are 
pushed to the real axis: (—00 + iO, — | + iO) |J [— f + z'0, a^jy] and (—00 — iO, — ^ — iO) |J [— | 
respectively. On [— ^,a 4 jv] and [— | — 2O, a 4 Ar + i] , w^?^ 7^ by continuity on a com- 
pact set. Finally, for all z G (— 00 + iO, — | + iO) , ^sh'(z,n) = — | and R(z,a) G R. So 
r\z%(1)) ^ f° r an 2 e (— 00 + 2O, — ^ + z'0) . The interval (— 00 — z'0, — ^ — z'0) is done 
similarly. So 7^ for all z G 7 00 , for any in a small neighborhood of yUo- 

□ 

Definition 3.13. 

We define the (finite) genus G = of the asymptotic solution of the semiclassical one 

dimensional focusing NLS with the initial condition defined through f(z, fi), as (finite) N G N 
if there exists an asymptotic solution of NLS through the solution h(z, fi) of RHP('~f oc , a, fi, /) 
with a = (a , cti, C(an+i)> such that h j^^ 7^ for all z G 7°°, and the signs conditions of 
h on 7 00 are satisfied: h G S , C(7°°). 
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Remark 3.14. The definition of the genus of the asymptotic solution can alternatively be 
seen as the genus of the model (limiting) Riemann surface. 

Theorem 3.15. (Preservation of genus) 

Suppose for fiQ, the genus of the asymptotic solution of NLS with initial condition defined 



through f(z,fio) in (15), is G(/U 



Then there is an open neighborhood of fio such that, for all fi in the neighborhood of fiQ, 
the genus of the asymptotic solution of NLS with initial condition defined through f(z,/i), is 
preserved G(fi) = G(fi ). 



Proof. Follows from Theorem 3.12 and Definition 3.13 



□ 

Corollary 3.16. Fix x andt > to, where to(x) is the time of the first break in the asymptotic 
solution. Then in some open neighborhood of fi — 2 the genus of the solution is 2. 

Proof. For fi = 2 and t > to(x) the genus is 2 for all x [2 



. By the Preservation of genus 
theorem |3.15[ the genus is preserved in some open neighborhood of fi = 2, including some 
open interval for fi < 2. □ 



3.4 Computations 



Figure |3| demonstrates comparison solutions of (70) and (63). The solutions are practically 
indistinguishable on the figure with the absolute difference less than 10 -3 for fi G [1,4], 
which includes a critical value (i = 2, the transition between (solitonless) pure radiation case 
(fi > 2) and the region with solitons (0 < fi < 2). 



4 Appendix 

4.1 Genus region 

It was shown in [23] that for all fi > and for all x, there is a breaking curve t = t (x) in 
the (x, t) plane. The region < t < t (x) has genus in the sense of genus of the underlying 
Riemann surface for the square root 



R(z) = \/{z- a Q )(z - a?T), 

where the branchcut is chosen along the main arc connecting «o and ot\ = ao, and the branch 
is fixed by R(z) — > —z as z — > +oo. The asymptotic solution of NLS is expressed in terms 
of «o = cto(x, t). 

All expressions in the genus region have simpler form. In particular: 

i / \ R{z,a ) r f((,fi)d( 

h{z,a ,fi) = —— <p — — — -, 58 

K(z,a ,fj.) = — <p — r, 59 

2tii h (C - z)R{(,a ) 
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x = 0.4 , t = 0.5 , n s [ 1 , 3 ] 






Re( a ) 



./.HO 




Figure 3: Comparison of /i evolution of a = (ao,a2,«4) using (70) and (63) 
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Oir 




0L\ 



«4 



Ot 2 




Q?5 0:3 ^1 



7m,0 



Figure 4: The jump contour in the case of genus and genus 2 with complex conjugate 
symmetry in the notation of [23] . 



and with a slight abuse of notation 

K{a Q ,n) := K(ao,a Q ,fi) 

and 



/(C,/*)dC 



27ri L (C - a )i?(C, Qf ) 



1 



dK OK 
djx dfx 2m L ((, - a )R(C, a ) 



(60) 



(61) 



Theorem 4.1. ([i-perturbation in genus 0) 

Consider a simple contour 70 = 7(^0) = [«o(/^o), «o(/^o)] an oriented simple arc with the 
distinct end points (a 7^ a ) and depending on parameter fx (see Figure^). Assume a (fx ) 
satisfies the equation 

K(a ,no) = 0, 



and f is given by (15). Let 7 = 7 (a, /i) be the contour of a RH problem which seeks a 
function h(z) which satisfies the following conditions 



h+(z) + h-(z) = 0, on 7, 
h(z) + f(z) is analytic in C\y, 



(62) 



Assume that there is a function h(z, do, /io) which satisfies (62) and suppose 
for all z on 7. 

Then the solution ao(fi) of the equation 

K (a ,fj) = 



h'(z,d ,n ) 
R(z,a ) 



7^0 



(63) 



and h(z,a(fi), fi) which solves (62) are uniquely defined and smooth in // in some neighbor- 
hood of /Jo ■ 

Furthermore 



dac 
dfj, 



(m) 



27ri^(a ,Ai) 



f'K) 



7(M) (C-aoM)-R(C>«o(A*)) 



(64) 



dh 

where z is inside 0/7. 



R(z,a Q ) 
2-ki 



7( M ) (C - 2)fl(C,ao(^)) 



(65) 



4.2 Genus 2 region 

In the genus 2 region (N = 1), with underlying Riemann surface for the square root 



R(z) = \/{z- a )(z - ai)(z - a 2 )(z - a 3 )(z - a A )(z - a 5 ), 

where the branchcut is chosen along the main arcs connecting o and a%, a 2 and o 4 , a 5 and 
a 3 ; and the branch is fixed by R(z) —> —z 3 as z — > +00. 
Taking into account the complex conjugate symmetry 



oil = oiq, a 3 



h{z) 



R(z) 



-d( + 



a 2 , a 5 
W 



«4. 



-d( + 



n 



-de 



(66) 
(67) 



2ttz LA (C " z)R(C) h m (C - z)R(C) K (C - z)R(C) 
where 7 m is a loop around the main arc [a 2 , an] IJ[ a 5; c^L an d 7c is a loop around the 
complementary arc [0:0,0:2] Ut^^i] ( see Fig- [p. 

Theorem 4.2. (^-perturbation in genus 2) 

Consider a simple contour 70 = 7(/xo) consisting of a finite union of oriented simple arcs 
7o — (U7m,j) U (UTc,j) ^i/i ^ e distinct arcs end points do = (00,^2,04) and depending on 
parameter /x (see Figure^. Assume a satisfies a system of equations 

K (a ,d Q ,iJ, ) = 0, 
K(a 2 ,a ,fi ) = 0, 
K (a 4 ,a ,/xo) = 0, 



and / is given by (15). Let 7 = 7(0, /x) fre i/ie contour of a RH problem which seeks a 



function h(z) which satisfies the following conditions 

h+(z) + h-(z) = 0, on 7 m>0 , 
/z + (z) + = 2W 7 , on 7 m> i, 

— /i-(^) = 2fi, on 7 Cj i, 
/i(z) + /(z) is analytic in C\y, 



(68) 



where Q = fl(a, /x) and = VT(a, /x) are rea/ constants whose numerical values will be 
determined from the RH conditions. Assume that there is a function h(z, do, /Xo) which 
satisfies (68) and suppose 7^ /or all z on 7. 

Taen the solution a = <3(/x) 0/ t/ie system 



K(ao,a Q ,fi) = 0, 
K(a 2 ,d ,fi) = 0, 
X (a 4 , d?o, /x) = 0, 



(69) 
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and h(z , a(n) , n) which solves (68) are uniquely defined and smooth in fi in some neighbor- 
hood of yUo ■ 

Moreover, fi(/i) = ana " W(aO = WO^CaOj A*) are defined and smooth in fi in 

some neighborhood of ji^. 
Furthermore 



da 3 

dh 
<9/i 



(z, a, x, t, ji) 



2-ni 



R(z, a 



2m 



|(c,aO 



7(m) 



(C-^)i2(C,a; 



where z is inside 0/7, 



dW 



x,t,n) ■■ 



1 

D 
1 





i 4£ 


J7m R(C) 


hm R(C) 






r d( 


£ 

K r(o 


K R(0 



where aj = aj(x,t,fi), R(() = R((, a(x,t, fi)), /(C) = f((,x,t,fi), f^Q = %(C,x,t, //) 



D = D(x,t,(j) 



R C) J* 



> R(0 Jj m R(0 
- _d(_ r (d(_ 

7c R(0 f% R(0 



(70) 



(71) 



(72) 



(73) 



and 



(74) 



References 

[1] Baik, J., Deift, P., Johansson, K., On the distribution of the length of the longest in- 
creasing subsequence of random permutations, J. Amer. Math. Soc. 12, no. 4 (1999), 
1119-1178. 

[2] Ceniceros, H., Tian, F.-R., A numerical Study of the semi- classical limit of the focusing 
nonlinear Schrddinger equation, Phys. Lett. A 306 (2002), no. 1, 25-34. 

[3] Belov, S., Tovbis, A., Venakides, S., Singular obstruction in the g-function mechanism, 
preprint. 

[4] Boutet de Monvel, A., Kotlyarov, V., Shepelsky, D., Focusing NLS equation: long-time 
dynamics of step-like initial data, Int. Math. Res. Not. IMRN 2011, no. 7, 16131653. 

[5] Buckingham, R., Venakides, S., Long-time asymptotics of the nonlinear Schrddinger equa- 
tion shock problem, Comm. Pure Appl. Math. 60 (2007), no. 9, 13491414. 

[6] Deift, P., Orthogonal polynomials and Random Matrices: a Riemann-Hilbert approach, 
Courant Lecture Notes in Mathematics, New York University, New York (1999). 



20 



[7] Deift, P., Kriecherbauer, T., McLaughlin, K. T.-R., Venakides, S., Zhou, X., Strong 
asymptotics of orthogonal polynomials with respect to exponential weights, Comm. Pure 
Appl. Math. 52 (1999), no. 12, 1491-1552. 

[8] Deift, P., Kriecherbauer, T., McLaughlin, K. T.-R., Venakides, S., Zhou, X., Uniform 
asymptotics for polynomials orthogonal with respect to varying exponential weights and 
applications to universality questions in random matrix theory, Comm. Pure Appl. Math. 
52 (1999), no. 11, 1335-1425. 

[9] Deift, P., Venakides, S., Zhou, X., New results in small dispersion KdV by an extension of 
the steepest descent method for Riemann-Hilbert problems, Internat. Math. Res. Notices 
6 (1997), 286-299. 

[10] Deift, P., Zhou, X., A steepest descent method for oscillatory Riemann - Hilbert prob- 
lems. Asymptotics for the mKdV equation, Ann. of Math. 137 (1993), 295-370. 

[11] Deift, P., Zhou, X., Asymptotics for the Painleve II equation, Comm. Pure and Appl. 
Math. 48 (1995), 277-337. 

[12] Duits, Maurice, Kuijlaars, Arno B. J., Universality in the two-matrix model: a Riemann- 
Hilbert steepest-descent analysis, Comm. Pure Appl. Math. 62 (2009), no. 8, 1076-1153. 

[13] Ercolani, N.M., McLaughlin, K.D.T-R., Asymptotics of the partition function for ran- 
dom matrices via Riemann-Hilbert techniques, and applications to graphical enumeration, 
Internat. Math. Research Notices 14, 755-820 (2003). 

[14] Fokas, A. S., Its, A. R., Kitaev, A. V., Discrete Painleve equations and their appearance 
in quantum gravity, Comm. Math. Phys. 142 (1991), no. 2, 313-344. 

[15] Fokas, A. S., Its, A. R., Kitaev, A. V., The isomonodromy approach to matrix models 
in 2D quantum gravity, Comm. Math. Phys. 147 (1992), no. 2, 395-430. 

[16] Lyng, C, Miller, P.D., The N-Soliton of the Focusing Nonlinear Schrddinger Equation 
for N Large, Comm. Pure Appl. Math. 60 (2007), no. 7, 951-1026. 

[17] Kamvissis, S.; McLaughlin, K. D. T.-R.; Miller, P. D., Semiclassical soliton ensembles 
for the focusing nonlinear Schrddinger equation, Annals of Mathematics Studies, 154. 
Princeton University Press, Princeton, N.J., 2003. 

[18] Miller P.D., Kamvissis, S., On the semiclassical limit of the focusing nonlinear 
Schrddinger equation, Phys. Lett. A 247 (1998), no. 1-2, 75-86. 

[19] Sabat, A. B., One- dimensional perturbations of a differential operator, and the inverse 
scattering problem, Problems in mechanics and mathematical physics, Nauka, Moscow, 
(1976), 279296. 

[20] Simon, B., Szego's theorem and its descendants, Princeton University Press, Princeton 
(2011) 



21 



[21] Tovbis, A., Venakides, S., Determinant Form of the Complex Phase Function of the 
Steepest Descent Analysis of RiemannHilbert Problems and Its Application to the Focus- 
ing Nonlinear Schrodinger Equation, IMRN no 11 (2009), 2056-2080. 

[22] Tovbis, A., Venakides, S., Nonlinear steepest descent asymptotics for semiclassical limit 
of integrable systems: continuation in the parameter space, Comm. Math. Phys. 295 
(2010), no. 1, 139160. 

[23] Tovbis, A., Venakides, S., Zhou, X., On semiclassical (zero dispersion limit) solutions 
of the focusing Nonlinear Schrodinger equation, Communications on Pure and Applied 
Mathematics, Vol. LVII, (2004) pp. 0877-0985. 

[24] Tovbis, A., Venakides, S., Zhou, X., On the long-time limit of semiclassical (zero dis- 
persion limit) solutions of the focusing nonlinear Schrodinger equation: pure radiation 
case, Comm. Pure Appl. Math. 59 (2006), no. 10, 1379-1432. 

[25] Tovbis, A., Venakides, S., Zhou, X., Semiclassical focusing nonlinear Schrodinger equa- 
tion I: inverse scattering map and its evolution for radiative initial data, Int. Math. Res. 
Not. IMRN 2007, no. 22, Art. ID rnm094. 



22 



